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ABSTRACT 

We present atmospheric parameters for 51 nearby F and G dwarf and subgiant stars uniformly dis¬ 
tributed over the —2.60 < [Fe/H] < +0.20 metallicity range that is suitable for the Galactic chemical 
evolution research. Lines of iron in the two ionization stages, Fe I and Fe II, were used to derive a 
homogeneous set of effective temperatures, surface gravities, iron abundances, and microturbulence 
velocities. Our spectroscopic analyses took advantage of employing high-resolution {R > 60 000) 
Shane/Hamilton and CFHT/ESPaDOnS observed spectra and non-local thermodynamic equilibrium 
(NLTE) line formation for Fe I and Fe II in the classical ID model atmospheres. The spectroscopic 
method was tested in advance with the 20 benchmark stars, for which there are multiple measure¬ 
ments of the infrared flux method (IRFM) effective temperature and their Hipparcos parallax error 
is less than 10%. We found NLTE abundances from lines of Fe I and Fe II to be consistent within 
0.06 dex for every benchmark star, when applying a scaling factor of 5'h = 0.5 to the Drawinian rates 


2 


Sitnova et al. 


of inelastic Fe+H collisions. The obtained atmospheric parameters were checked for each program 
star by comparing its position in the log g - Tl-fF plane wi th the theoretical evolutionary track of 
given metallicity and a-enhancement in the lYi et al.l (l2004l g rid. Our final effective temp eratures 
lie exactly in between the Tirfm scales of lAlonso et al.l ([l996af) and ICasagrande et al.l (|2011f) . with a 
mean difference of +46 K and —51 K, respectively. NLTE leads to higher surface gravity compared 
with that for LTE. The shift in log g is smaller than 0.1 dex for stars with either [Fe/H] > —0.75, or 
TeB < 5750 K, or log g > 4.20. NLTE analysis is crucial for the VMP turn-off and subgiant stars, 
for which the shift in log g between NLTE and LTE can be up to 0.5 dex. The obtained accurate 
atmospheric parameters will be used in the forthcoming papers to determine NLTE abundances of 
important astrophysical elements from lithium to europium and to improve observational constraints 
on the chemo-dynamical models of the Galaxy evolution. 

Subject headings: Stars: abundances - Stars: atmospheres - Stars: fundamental parameters 


1. INTRODUCTION 

In recent decades the construction of large telescopes, 
the developement of efficient spectrometers, and progress 
in stellar atmosphere and synthetic spectrum numerical 
modelling provided a considerable improvement in 
the quality and quantity of stellar abundance deter¬ 
minations. Elemental abundances for the FGK-type 
stars provide important clues to understand the main 
processes at play in formation and evolution of the Milky 
Way (MW). Studies of stellar samples with metallicity 
[Fe/H]^ > —1 showed how trends in various chemical 
elements can be used to learn t he chemodynami- 
cal e v olution of the Ga lactic disk (lEdvardsson et ahl 
1199,11: iGhen et ^ l2000h and to r esolve the thick 
disk and thin disk (iGratto n et al.l 11996: iFuhrmann 
19981 l2004 [Mashonkina &: Gehrenl 120001: Reddv et al. 


iMishenina et al.l l2004t I Zhang fc Zhaol l2006t 


Adibekvan et al.l 1201, It iBensbv et al.l I20l4 and refer¬ 


ences therein). Studies of very metal-poor (VMP, [Fe/H] 
< —2) stars unraveled the maj or processes of chemical 
enrichment of the Milky Way (E cWilliam et al.1 1199, 51 
Gavrel et al.ll2004 IZhang fc Zhaol 1200, 4 iBonifacio et al l 


200911 . 


Elemental abundances of stars with different metallic- 
ities are the mai n constraint for Galacti c chemodynami¬ 
cal models ( e.g. iGhiappini et al.1 (|2001lj ; iRomano et al.1 
(|2010lj ; iKobavashi et al.l ( 2011h j. For the Galactic chem¬ 
ical evolution research it would be useful to deal with a 
homogeneous set of accurate stellar abundances in wide 
metallicity range, from supersolar down to extremely 
low iron abundances. Large high-resolution spectro¬ 
scopic surveys were proposed to increase the statistics 
of observations and to improve a homogeneity of de¬ 
rived stellar chemical abundances over wide metallicity 
range. The Apache Point Obse rvatory Galactic Evo - 
lution Experiment (APOGEE) (jZasowski et al.l 1201,111 . 
with its 10® red giants observed in the near-infrared 
H-band with a spectral resolving power of i? ~ 22 500, 
enables users to address numerous Galactic structure 
and stellar populations issues. The Gaia-ESO Survey 
consortium (iGilmore et ahl l2012t iRandich et al.l 120131 : 
LSmilianic et al.ll2014lj is obtaining high-quality spectro¬ 
scopic data for about 10® stars using FLAMES at the 
VLT. Spectra for a million stars will be acquired by the 
coming Galactic Arch aeolog y with HERMES (GALAH) 
survey (iDe Silva et al.l[20T^ . 

Electronic address: sitnova@inasan.ru 

^ In the classical notation, where [X/H] = log(A^x/-^H)star — 
log(A^x/A^H) Sun • 


We initiate a new project of derivjng a homogeneous 
set of stellar atmosphere parameters and chemical abun¬ 
dances for the Galactic field FGK-type stars in the metal- 
licity range —3 < [Fe/H] < +0.3 that is suitable for a 
systematic research of Galactic chemical evolution. By 
employing high-resolution spectral observations, deriv¬ 
ing accurate stellar atmosphere parameters, and treating 
the non-local thermodynamic equilibrium (NLTE) line 
formation for the key chemical species, we attempt to 
push the accuracy of the abundance analysis to the points 
where the trends with metallicity could be meaningfully 
discussed. The selected stars are uniformly distributed in 
metallicity and they can serve as a calibration sample for 
the existing and coming large-scale stellar surveys, such 
as the LAMOST Experim ent for Galactic U nderstanding 
and Exploration ILEGUE lDeng et'nil2012i) . designed to 
determine from low-resolution spectra elemental abun¬ 
dances for hundreds of thousands to millions of stars, 
spread over much larger distances than ever before. 

The goal of the present paper is the determination of 
precise atmospheric parameters, i.e. the effective tem¬ 
perature, TeB, the surface gravity, log g, the iron abun¬ 
dance, [Fe/H], and the microturbulence velocity, ^t, for 
the selected sample of, presumably, dwarf stars. The 
next papers in the series will concern with calculations 
of the NLTE abundances for many chemical elements 
from Li to Eu and analysis of the obtained abundance 
trends. Taking advantage of the Shane/Hamilton and 
GFHT/ESPaDOnS high-resolution observational mate¬ 
rial with sufficient spectral coverage and the NLTE line 
formation for Fe I-Fe H, we determined stellar atmo¬ 
sphere parameters spectroscopically from lines of iron in 
the two ionization stages, Fe I and Fe H. 

The spectroscopic methods are, in particular, useful 
for evaluating Tqb and log g of the VMP stars that are 
mostly distant, and, therefore, their temperatures can¬ 
not be reliably derived from photometric colours due to 
the uncertainty in the interstellar reddening data and, 
at present, surface gravities cannot be calculated from 
the trigonometric parallax because it is either rather un¬ 
certain or non-measurable. We look forward to seeing 
soon accurate parallaxes from Gaia^. In the literature 
there is no consensus on a validity of the excitation tem¬ 
peratures, Texc, derived through forcing no dependence 
of the Fe I line abundance on the excitation energy of 
the lo wer level, Ee ^ e- Fo r their sample of VMP cool 
giants iGavrel et al.1 (|2004 j reached in LTE good agree¬ 
ment between the photometric and spectroscopic tem- 

^ http://sci.6sa.int/gaia/ 
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peratures after they excluded strong lines with i?exc = 
0 eV. A similar approach (LTE, but not cutting Eexc = 
0 eV) yielded the same outcome for t he VMP near-main- 
sequence stars in iFrebel et al.l (|2013ll . while up to several 
hundred degrees cooler excitation temperatures than the 
p hotometric ones wer e found for the cool giants. 

iRuchti et al.l (|2013ll determined surface gravities of the 
stellar sample selected from the RAVE surv ey, using a 
grid o f the NLTE abundance corrections from lLind et al.l 
(|2012ll . For common in the NLTE and LTE analysis, 
they found systematic biases in log g of up to 0.2 dex in 
the metallicity range — 2 < [Fe/H] < —0.5 and up to 
0.3 dex for the more metal-poor stars. 

The NLTE effects for lines of Fe I were accounted for 
by iBensbv et al.l (|2014f) in their study of 714 F and G 
dwarf stars, by applying the NLTE abundance correc¬ 
tions interpolated in the grid of iLind et al.l (l2012ll . They 
found minor shifts between NLTE and LTE, with ATeff = 
-12 ±28 K, Alog g = ±0.012 ± 0.059, and Aloge(Fej = 
—0.013 ± 0.016, on average. This is because the major¬ 
ity of their stars are close-to-solar metallicity and mildly 
metal-deficient ones, with [Fe/H] > —1.2, and the NLTE 
abundance corrections are sma ll, at the level of few units 
of h undredth, a s pred icted by iMashonkina et all (|201ir) 
and iLind et ^ (l2012fl for this metal licity range , fur¬ 
thermore, the NLTE calculations of iLind et ^ ()2ni2f l 
resulted in small departures from LTE for Fe I until 
the extremely low metallicities because they were per¬ 
formed assuming a high efficiency of the Fe±H coll i- 
sions and using the Drawinian rates (|Drawinlll968lll96^ . 
For nearby stars with very go o d Hip parcos parallaxes 
(Att/tt < 0.05), IBensbv et ahl (|2014ll found that essen¬ 
tially all stars with log g > 4.2 and < 5650 K do 
not show ionization equilibrium between Fe I and Fe II, 
when deter mining the su r face g ravity from Hipparcos 
parallaxes. IBensbv et al.1 (|2014f l suggested that classical 
plane-parallel (ID) models have limitations and cannot 
properly handle excitation balance and/or ionization bal¬ 
ance, and they applied, therefore, empirical corrections 
to the atmospheric parameters from ionization balance. 
In this study we check with our stellar sample, whether 
spectroscopic methods of and log g determination 
from lines of Fe I and Fe II have any limitations and 
what they are. 

The paper is structured as follows. In Sect. [5] we de¬ 
scribe the stellar sample, observations, and their reduc¬ 
tion. Kinematical properties of the selected stars are cal¬ 
culated in Sect. |3l They are used to identify a member¬ 
ship of individual stars to the galactic stellar populations. 
Section [4] concerns with the methodical issues. Stellar ef¬ 
fective temperatures, surface gravities, iron abundances, 
microturbulence velocities, and masses are derived in 
Sect, m The obtained results are discussed in Sect. 1 
Our conclusions are presented in Sect. 0 

2. STELLAR SAMPLE, OBSERVATIONS, SPECTRA 
REDUCTION 

The stars were select e d from the [Fe/H] catalogue of 
ICavrel de Strobel et al.l (|2001ll based on the following 
criteria. 

1. The stars have a declination of J > —20° to be 
observed at the nothern sky. 

2. The selected stars should cover as broad as possi¬ 


ble metallicity range and be uniformly distributed, 
with 2-3 stars in each 0.1 dex metallicity interval. 

3. The stellar sample should be homogeneous in tem¬ 
perature and luminosity. We selected the F-G-K 
dwarfs and subgiants based on the literature data. 

4. Binaries, variables, and stars with any chemical pe¬ 
culiarity (carbon-enhanced stars, low [a/Fe] stars, 
etc.) were excluded. 

5. For testing purposes the stellar sample should in¬ 
clude a dozen of well-studied stars of various metal¬ 
licities, for which their stellar atmosphere parame¬ 
ters Teff and log g were determined in the literature, 
presumably, from the non-spectroscopic methods. 

Thus, 50 stars, in total, were selected, and they cover 
the —2.6 < [Fe/H] < ±0.3 metallicity range. 

Spectra of 48 stars were obtained using the Hamilton 
Echelle Spectrograph mounted on the Shane 3-m tele¬ 
scope of the Lick observatory during the two observa¬ 
tion runs March 15-16, 2011 and January 5-11, 2012. 

Most stars were observed, at least, twice. The resolving 
power IS R — A/AA ~ 60 000, and the spectral coverage 
is 3700 A to 9300 A. The signal-to-noise ratio {S/N) at 
5500 A is higher than 100 for most stars. The excep¬ 
tion is BD -04° 3208, where S/N{X = 5500A) ~ 70. It 
is worth noting that the spectra suffer from the fringing 
effect in the infrared band. 

For two stars, we used the spectra observed with 
the ESPaDOnS echelle spectrograph at the 3.6-m tele¬ 
scope of the Canada-France-Hawaii (CFH) observatory 
in queued service observing mode during several nights 
in 2011 and 2012. The spectra cover 3700 A to 10450 A. 

They have i? ~ 81 000 and S/N ~ 150 at 5030 A and 
higher value of about 200 at 8090 A. 

The Shane/Hamilton spectra were reduced with the 
IRAF package following standard procedures that in¬ 
clude the aperture definition, flat-field correction, back¬ 
ground subtraction, one-dimensional spectra extraction, 
wavelength calibration, and continuum normalization. 

The wavelength calibration was performed using the 
lamp with the titanium cathode in argon environment. 

Such a lamp is rare in the astronomical observations, and 
an appropriate list of t he Ti and Ar lines was compiled 
bv iPakhomov fc Zhaol ([20131 1. The uncertainty in wave- 
length calibration is estimated to be 0.006 A. 

The CFHT/ESPaDOnS spectra were reduced with 
the special reduction package Libre-ESpRIT^, which in¬ 
cludes two steps. The first one performs a geometri¬ 
cal analysis from a sequence of calibration exposures, 
and the second step achieves spectrum optimal extrac¬ 
tion in itself, using the geometrical information derived 
in the first step. The reduction procedure returns the 
normalised spectra. 

To improve the statistics of the VMP stars, we in¬ 
clude the well-studied halo star HD 140283 in our sam¬ 
ple. Its high-quality observed spectrum was taken from 
the ESO UVESPOP survey (iBag-nulo et al.ll200^ . For 
seven stars our observational material was complemented 
with the data from different sources. We employed the 

® http://www.cfht.hawaii.edu/Instruments/Spectroscopy/Espadons 
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Fig. 1.— Toomre diagram for the investigated stars from the 
perspective of the local standard of rest. Dashed curves delineate 
constant peculiar space velocities Upec = {U^ + = 85, 

180, and 400kms~^. Different symbols show the thin disk (open 
circles), thick disk (filled circles), and halo (open rhombi) stars. 

spectra obtained by Klaus Fuhrmann with the fiber op¬ 
tics Cassegrain echelle spectrograph FOCES at the 2.2-m 
telescope of the Calar A lto Observatory in 1996 to 2000 
(|Fuhrmannlll998L I2004D . The star BD—04° 3208 was ob¬ 
served using UVES/VLT in April 2001 within our project 
67.D-0086A (PI: T. Ge hren). Details of sp ectra reduc¬ 
tion were described by iMashonkina et aTT (l2003l) . For 
HD 49933 its high-quality spectrum was observed with 
the HARPS spectrograph at the 3.6-m ESO La Silla tele- 
scope. Details of s pectr a reduction were described by 
iR.vabchikova et aO ()2009[ ) . 

This paper deals with 51 stars, in total. The investi¬ 
gated stars together with characteristics of the observed 
spectra are listed in Table[2] 

3. KINEMATIC PROPERTIES OF THE SELECTED STARS 

The stellar kinematics is closely related to stellar pop¬ 
ulations in the Galaxy, and it is commonly applied to 
identify a membership of given star to the galactic stel¬ 
lar populations, namely the thin disk, the thick disk, and 
the halo. 

The galactic space velocity components U, V, W 
were calculated using the e quations and formalism of 
[Johnson fc SoderblomI (119871) . They were defined with 
respect to the local standard of rest (LSR), adopting 
the standard s olar motion (U, W) = (10.00, 5.25, 
7.17) kms“^ of iDehnen fc BinnevI (I1998D . In computa¬ 
tions of the X, Y, Z-coordinates, we used the best cur¬ 
rent estimate of the Galactocentric distance of the Sun 
Rg = 8.0 kpc, which was inferr ed from a com parison of 
different statistical techniques (iMalkinI 120131) . The par¬ 
allaxes and proper motions were take n from the updated 
version of the Hipparcos catalogue (Ivan Leeuwenl[200^ 
and the radial velocities f rom the Hipparcos Input cat¬ 
alogue (jTuron et al.lll993l) . The obtained results are pre¬ 
sented in TablelH The star BD—13° 3442 is missing in 
the Hipparcos catalogue, and it was assigned to the 
galactic halo based on its very low Fe abundance, with 
[Fe/H] = —2.62 (Sect. [5]), and high proper motion as in¬ 
dicated by t he Simbad databa se^. 

Foll owing iFuhrm annI ( 2000l) and later stu dies (for ex¬ 
ample [ChS(^^ni2Q3 Bensbv et al.ll20l6l) . we identify 

http://simbad.u-strasbg.fr/simbad/ 


the stars with peculiar space velocities with respect to 
the LSR in the range 85 kms“^ < Vpec < 180 kms“^ as 
belonging to the thick disk stellar population, the stars 
with Vpec < 85 kms“^ to the thin disk, and all the halo 
stars of our sample have Vpec > 200 kms“^ (Fig.[I|). The 
star G090-003 is not shown in Fig.[Tl because its pecu¬ 
liar space velocity, Vpec = 847 ± 1580 kms“^, is rather 
uncertain due to large error of the Hipparcos paral¬ 
lax, TT = 1.12 ± 2.14 mas. For six stars with peculiar 
space velocities close to Vpec = 85 kms“^, their member¬ 
ship to either the thin disk or the thick disk cannot be 
decided us i ng on ly the kinematic properties. Following 
iFuhrmaiml (119981 ). we also took into account their chem¬ 
ical properties, as determined in Sect.jSJ For example, 
HD30562, with Vpec = 93 kms“^, was assigned to the 
thin disk because it has a supersolar Fe abundance, with 
[Fe/H] = 0.17, and it does not reveal Mg enhancement 
relative to Fe. 

4. METHOD OF CALCULATIONS 

Stellar parameters of the selected sample were deter¬ 
mined in this study applying different methods, includ¬ 
ing the spectroscopic one based on the NLTE analysis of 
lines of Fe I and Fe H. This section describes calculations 
of the theoretical spectra. 

4.1. Codes and model atmospheres 

The statistical equilibrium (SE) of Fe I-Fe H was 
cal culated with a comprehe nsive model atom treated 
by IMashonkina et ahl (120111) . The main source of the 
uncertainty in the NLTE results for Fe I is poorly 
known inelastic collisions with hydrogen atom s. We 
employed the formula of iDrawi ^ IT^ , as im- 

plemented bv iSteenbock fc Holwegeil (119841) . for allowed 
b — b and b — f transitions and a simple relation be¬ 
tween hydrogen and electron collisional rates, Ch = 
Ce\ /(me/rnH)NH/ Ne . for forbidden transitions , follow¬ 
ing iTakedal (jl994l ). IMashonkina et all (|20111 ) recom¬ 
mended to scale the Drawin rates by a factor of S'h = 
0.1 based on their analysis of element abundances de¬ 
rived from the two ioniza tion stages, Fe I and Fe H, in 
the selected VMP stars. iBergemann et al.l (120121) esti¬ 
mated a larger value of S'h = 1. In this study, was 
further constrained using the benchmark stars. 

The coupled radiative transfer and SE equations were 
solved with a revised v ersion of the DETAIL code 
(iButler fc Giddin^ll98^ . The update was presented by 
IMashonkina et all (|20111l . The obtained dep arture coeffi¬ 
cient s were then used by the codes binmag3 (jKochukhovI 
IMol) . synthV-NLTE (Vadim T symbal, p rivate commu- 
nicatiqn, bas ed on SynthV code lTsvmballll99CT and SIU 
(|Reetdll9911) to calculate the synthetic line profiles. 

We used t he _ MARGS model structures 

(|Gustafsson et al.l 120081) ^. which were interpolated 
for given Teg, log g, and [Fe/H] using a FORTRAN- 
based routine written by Thomas Masseron®. 

4.2. Selection of iron lines 

The investigated lines were selected according to the 
following criteria. 

® http://marcs.astro.uu.se 

® http://marcs.astro.uu.se/software.php 
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(i) The lines should be almost free of visible/known 
blends in the solar spectrum. 

(ii) For Fe I the linelist should cover a range of exci¬ 
tation energies of the lower level, Fexc, as large as pos¬ 
sible to investigate the excitation equilibrium of neutral 
iron. On the other hand, the low-excitation (Fexc < 
1.2 eV) lines in metal-poor stars often appear to have 
higher abundances t han lines with higher excitation en - 
ergy (see, for examole fCavrel et al.ll2004lLai et al.ll2008ll . 
The calculations with ab initio 3D, time-dependent, hy- 
drodynamical model atmospheres of cool stars predict 
that the 3 D effects for lines of Fe I are strongly Fpx r 
dependent (|Collet et al.ll2007l : iDobrovolskas et al.feoi^ . 
For example, the 3D abundance corrections amount to 
—0.8 dex and —0.25 dex for the Eg^c = 0 lines of Fe I 
in the models Tgg / log g / [M/H] = 5860/4/—2 and 
5850/4/—1, respectively, while they are —0.2 dex and 
—0.07 dex for the Eg^g = 2 eV lines and close to 0 for 
the Fexc = 4 eV lines (|Dobrovolskas et 811120131 1. To 
minimise a possible influence of the 3D effects on Tgg 
determinations, we did not consider lines of Fe I with 
Fexc < 2 eV. 

(iii) For both close-to-solar metallicity and VMP stars, 
lines of various strength have to be present in the spec¬ 
trum to evaluate stellar microturbulence velocity 

For individual stars, we also avoided using of saturated 
lines with the equivalent width EW > 180 mA to min¬ 
imise an influence of possible uncertainty in the van der 
Waals damping constants on final abundances. 

The selected lines are listed in Tabled along with their 
atomic parameters. For neutral iron we employed exp er¬ 
imental (?/-val u es fro m iBlackwell et al.l (fl979L Il982al lbll ; 
lO’Brian et al.l (|1991ll . and iBard et al.l (1199111. Two 
sets o f osci llator strengths from iMelendez fc Barbuvl 
(j2009t) and iRaassen fc Uvlin^ (|1998ll were inspected 
for Fe II. Van der Waals broadening of the selected 
lines is accounte d for using th e mos t accurate data 


as provided by iBarklem et al. 

j200Cf) for Fe I and 

[Barklem & Asoelund-Johansson 

pio^ for Fe H. 


4.3. Sun as a reference star 

To minimise the effect of the uncertainty in 5 /-values 
on the final results, we applied a line-by-line differen¬ 
tial NLTF and LTF approach, in the sense that stellar 
line abundances were compared with individual abun¬ 
dances of their solar counterparts. The solar flux obser- 
vations were taken f rom the Kitt Peak Solar Flux Atlas 
(iKurucz et al.l 119831 . The calculations were performed 
with the MARCS model atmosphere 5777/4.44/0. A mi¬ 
croturbulence velocity of 0.9kms“^ was adopted. 

The solar LTF and NLTE abundances from individ¬ 
ual lines are presented in TableSl The NLTE calcula- 
toins were performed using = 0.5 (see sect. B- For 
the Fe I lines, the solar mean LTF and NLTE abun¬ 
dances are log Alte = —4.57 ±0.09 and log Anlte = 
—4.56 ± 0.09. Hereafter, Ax = Nx/Ntot is the ele¬ 
ment abundance taken relative to the total number of 
atoms. For lines of Fe II, the NLTE abundance cor¬ 
rections do not exceed O.OI dex, in absolute value, and 
the solar mean abundance amounts to log A = —4.56 ± 
0.05 and log A = —4.50 ± 0.05 , w hen using ff/-values 
from IMelendez fc Barbuvl (|2009ll and IRaassen fc Uvlingsl 
(|1998[) . respectively. Hereafter, the statistical abundance 
error is the dispersion in the single line measurements 


about the mean: a = •\/S(a: — Xi)^/{N — 1), where N is 
the total number of lines used, x is their mean abundance, 
Xi the individual abundance of each line. It is worth not¬ 
ing that an uncertainty of 0.06 dex in abundance from 
the Fe H lines leads to an 0.12 dex uncertainty in log g 
for solar-type stars. This is crucial for surface gravity de¬ 
terminations from the Fe I/Fe H ionization equilibrium. 

5. DETERMINATION OF STELLAR PARAMETERS 

To derive atmospheric parameters of our stellar sam¬ 
ple, the following strategy was applied. First, we selected 
the stars, for which their effective temperatures and sur¬ 
face gravities were determined reliably using the non- 
spectroscopic methods, namely Teg from the infrared flux 
method (IRFM) and log g from the well-known relation 
between log g, stellar mass M, Tgg, and absolute bolo- 
metric magnitude MboP 

[ 5 ] = [M] ± 4[Teff] ± 0.4(Mboi - Mboio). (1) 

Here, square brackets denote the logarithmic ratio 
with respect to the solar value. The Mboi magni¬ 
tudes were obtain ed using the V magnitudes from 
lOlsenI (I1983L 1199311 . th e revised Hipparcos parallaxes 
of Ivan LeeuwenI (120071). a.nd th e bolometric corrections 
('BCl from lAlonso et al.l (|1995ll . For the Sun the abso¬ 
lute visual magnitude My© = 4.87 and BC© = —0.12 
were adopted. The star’s mass was estimated from its 
position in the Mboi - log Tgg diagram by interp olating 
in the Dartmouth isochrones (iDotter et al.ll200^ . “Re¬ 
liably” means that multiple measurements of the IRFM 
temperature are available and the uncertainty in the HiP- 
PARCOS parallax based surface gravity, log gHip, is less 
than 0.1 dex. This implies a parallax error of smaller 
than 10%. The 20 stars in our sample meet these re¬ 
quirements. They are listed in Table[5]as the benchmark 
stars. 

Atmospheric parameters of the remaining stars were 
derived or improved spectroscopically from the NLTE 
analysis of lines of Fe I and Fe H. 

In the stellar parameter range, with which we concern, 
lines of Fe I are weakened towards higher Tgg, result¬ 
ing in higher derived element abundance, while they are 
nearly insensitive to variation in log g. In contrast, lines 
of Fe H are only weakly sensitive to variation in Tgg, 
while they are weakened with log g increasing. If Tgg is 
fixed, the surface gravity is obtained from the require¬ 
ment that abundances from lines of Fe I and Fe H must 
be equal. One needs to be cautious, when deriving both 
Tgg and log g from Fe I and Fe H. Any shift in Tgg 
leads to shifting log g in the sense: an increase in Tgg 
leads to an increase in log g. Therefore, in this study we 
used non-spectroscopic (IRFM) determinations of Tgg, 
where available, and for every star its final parameters, 
Teff/log g/[Fe/H], were checked with the corresponding 
evolutionary track. 

For each benchmark star, abundances from lines of Fe I 
and Fe H were derived under various line-formation as¬ 
sumptions, i.e., NLTE with = 0.1, 0.5, 1 and LTF, 
and we investigated which of them leads to consistent 
element abundances from both ionization stages. Our 
calculations showed that the departures from LTF are 
small for the [Fe/H] > — 1 stars. Indeed, in this metal- 
licity range the NLTE abundance correction, Anlte = 
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TABLE 1 

Abundance differences Fe I - Fe II in the 

SELECTED THREE STARS FOR DIFFERENT LINE 
FORMATION SCENARIOS. 


HD 

LTE 

Sh = 1 

NLTE 
Sh = 0.5 

II 

P 

84937 

-0.06 

-0.02 

0.00 

0.11 

94028 

-0.06 

-0.05 

-0.04 

-0.01 

140283 

-0.06 

-0.04 

-0.02 

0.10 


log^NLTE — log^LTE, does not exceed 0.04 dex for 
= 0.5, and the difference in NLTE abundances between 
= 0.1 and 0.5 is smaller than 0.04 dex. The three 
of six [Fe/H] < —1 stars are cool, with Tee < 5400 K, 
and they are not useful for constraining S'h because of 
small Anlte < 0-02 dex. Tabled] lists the abundance 
differences between Fe I and Fe II from the calculations 
with different S'h for the three benchmark stars with 
the largest NLTF effects. For every star, the accept¬ 
able abundance difference is achieved with Sh =0.5 and 
1. Hereafter, we choose Sh = 0.5 to perform the NLTF 
calculations for Fe I-Fell. 

Spectroscopic stellar parameters were obtained in the 
iterative procedure. For each star, our initial guess was 
Tee = Tirfm, if available, or Tee = Tb^y and log g = 
log gHip- Effective temperature and surface gravity were 
corrected until the ionisation equilibrium between Fe I 
and Fe II and the Fe I excitation equilibrium were ful¬ 
filled. 



[Fe/HJ 

Fig. 2.— Top panel: differences in th e IRFM effective 
tempe ratures of the investigated st ars be tween ICasagrande et ^ 
1120111 . Cll) and lAlonso et al.| l|1996a|. A96) (filled c ircle^ 

and between Cll and IGonzalez HerndrSez &: Bonifacio! 1^20091 . 
GB0 9) (open circles). M iddle panel: differences 'riRFM(A96) - 
Ti, llAlonso et al.lll996bh . Bottom panel: differences Tjrpm(C ll) 
- (ICasagrande et al.ll201CII . CIO). 


5.1. Effective temperatures 

Although there are multiple measurements 
of the IRFM effective temperatu re for in¬ 
dividu al stars of our sample lAlonso et ahl 

^96^ _ IGonzalez Hern andez fc Bonifaciol 120091 : 

Casagrande et al.l I 20 ldr ~ l 201 l[ ). we found no study 
that provides the data for every investigated star. 
Figured) (top panel) displays the differences between 
different sources for the stars in comm on with this work. 
It can be seen that the Tee scale of ICasagrande et ah 
(20141 Cll) is hotter than that of lAlonso et ah 
Jl996aL A96), with ATirfm(C11 - A96) = 114 ± 
62 K for 25 stars. The differe nc e betw een Cll and 
IGonzalez Hernandez fc Bonifaciol (|2009l GB09) is 

smaller, AriRFM(Cll - GB09) = 32 ± 83 K for 26 stars. 
Here, we did not count the two outliers HD 34411 and 
HD 142373, with large uncertainties in Tirfm of 479 
K and 342 K, respectively, that were caused by using 
saturated 2MASS photometry. 

We calculated four sets of the photometric temper¬ 
atures from the b — y and V — K colour indices us- 
ing two dif ferent calibrat i ons ofICasagrande et al.l (|2010L 
CIO) and lAlonso et al.l (|1996blL ' The V magnitudes 
and b — y colours were taken from lOlsenI (| 19831 Il993l) 
for 45 stars, in to tal, and the K magn itudes from the 
2MASS catalogue (ISkrutskie et al.l[2006D for 24 stars, in 
total. Figured) display the differences between Tirfm 
and Tb-y for both effective temperature scales. We 
found that Tb-y and Tv-k agree well with those from 
the IRFM, when using the calibrations of CIO, with 
Tirfm(CII) - Tb-yiClO) = 10 ± 46 K (42 stars) and 
Tt rfm(CII) - Tl f y(Cl O) = 22 ± 68 K (21 stars). For 
the lAlonso et al.l (|1996bl) calibration, the differences are 


larger, with Tirfm (A96) - Tb-y = 41 ± 78 K (25 stars) 
and Tirfm(A96) — Ty-j^ = 72 ± 58 K (12 stars), al¬ 
though the statistics is small for the V — K photometry. 

For the benchmark stars, we aimed to have a homoge¬ 
neous set of Teff based on a single source of the IRFM 
temperatures. The only study that provides Tirfm for 
every benchmark star is GB09. Using their Tirfm, we 
could not achieve the Fe I/Fe H ionization equilibrium 
and the Fe I excitation equilibrium for half benchmark 
stars, independent of applying either LTE or NLTE ap¬ 
proach. No preference was, therefore, given to any source 
of Tirfm- Instead, for each star its temperature from 
A96, GB09, CIO, and Cll, where available, and the cor¬ 
responding log gHip value were checked with the differ¬ 
ence in NLTE abundances between Fe I and Fe H and 
a slope of the [Fe/H] versus Eexc plot for Fe 1. For six 
stars, corrections up to 50 K were applied to the most ap- 
propiate Tirfm to obtain the final effective temperature 
(TableO. 

5.2. Surface gravities, metallicities, and 
microturhulence velocities 

For each benchmark star, the NLTE abundances from 
the two ionization stages, Fe I and Fe H, were found to 
be consistent within 0.06 dex. Hereafter, 0.06 dex was 
considered as an admissible difference between Fe I and 
Fe H, when deriving spectroscopic stellar parameters of 
the remaining, non-benchmark, stars. For given star, 
we started from checking the photometric temperature, 
Tirfm or Tb-y, and the Hipparcos parallax based grav¬ 
ity, log gHip- Effective temperature and surface gravity 
were allowed to vary within the corresponding error bars. 
The procedure was iterated until the NLTE abundance 
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Fig. 3. — Absolute (top panel) and differential (bottom panel) 
NLTE abundances from lines of Fe I (open circles) and Fe II (filled 
circles) in HD 22484. 
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Fig. 4.— Absolute (top panel) and differential (bottom panel) 
NLTE abundances from lines of Fe I (open circles) and Fe II (filled 
circles) in HD 84937. 


difference Fe I - Fe II is getting smaller than 0.06 dex, the 
excitation trend for Fe I disappears, and lines of different 
equivalent width give consistent iron abundances. 

We found that the differential approach largely re¬ 



Wavelength 


Fig. 5.— Theoretical NLTE profiles of Fe I 5383 A computed 
with logCe = —3 0.370 IIBarklem et al.ll2000l. dashed curve) and 
logCg = —31.095 l^mr^ d mon cont inuous curve) compared to 
the solar spectrum ([Kuru^et^l1ll984l . bold dots). Everywhere in 
the calculations, log Ape = —4.54. To fit the solar line profile with 
logCg = —30.370, one needs to reduce the iron abundance down 
to logApe = —4.72. 

moves a line-to-line scatter for the [Fe/H] > —1.5 stars. 
For example, when moving from the absolute to the dif¬ 
ferential NLTF abundances, the dispersion reduces from 
tr = 0.10 dex to tr = 0.04 dex for lines of Fe I and 
from cr = 0.08 dex to u = 0.05 dex for lines of Fe II in 
HD 22484 (6000/4.07/0.01, Fig.|31). For more metal-poor 
stars, the differential approach is less efficient in remov¬ 
ing a line-to-line scatter. For example, for HD 84937 
(6350/4.09/—2.12) the scatter of data reduces for Fe H, 
but not Fe I lines (Fig.|4|). When using the Fexc > 2 eV 
lines of Fe I, we obtained a = 0.035 dex and 0.078 dex 
for the absolute and differential abundances, respectively. 
This is, probably, due to the uncertainties in the van der 
Waals damping constant, Cq- The lines, which can be 
measured in the VMP stars, have strong van der Waals 
broadened wings in the solar spectrum (for example Fe I 
5383 A, Fig. [5]). For such a line, the uncertainty in the 
derived solar abundance is contributed from the uncer¬ 
tainties in both gf- and Ce-values. When dealing with 
the VMP stars, the differential analysis is able to cancel 
the uncertainty in gf, but not Cq. 

Although our final results are based on using the 
li'exc > 2 eV lines of Fe I, for most stars we also checked 
abundances determined from the lower excitation lines. 
We found no abundance difference between the low and 
high excitation lines for close-to-solar metallicity stars. 
For the most metal-poor stars, such as HD 84937, the 
absolute abundances reveal an excitation trend, but it 
disappears for the differential abundances (Fig.|4l). 

We did not determine stellar parameters under the 
LTE assumption and cannot evaluate the differences in 
log g between using LTE and NLTE. However, the NLTE 
effects can be estimated based on our LTE and NLTE cal¬ 
culations for Fe I-Fe H and the sensitivity of iron lines to 
variation in log g. Since the departures from LTE lead 
to weakened lines of Fe I, but they do not affect lines of 
Fe H until the extremely low metallicities, the ionization 
equilibrium Fe I/Fe H is achieved for higher gravities in 
NLTE than in LTE. The shifts in log g increase towards 
lower metallicity and surface gravity. For our program 
stars they may be up to 0.1 dex at [Fe/H] > —1.5, up 
to 0.2 dex for —2.2 < [Fe/H] < —1.8 and reach 0.45 dex 
for the most MP star of our sample BD —13° 3442, with 
[Fe/H] = -2.62. 

5.3. Checking atmospheric parameters with evolutionary 

tracks 
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Fig. 6.— Top panel: investigated thick disc (filled circles) and 
halo (diamonds) stars compared with the evolutionary tracks of M 
= 0.75 Mq and [Fe/H] from —2.75 to —0.75 (from left to right), 
with a step of 0.5 dex. Bottom panel: the thin disc stars (open 
circles) compared with the evolutionary tracks of the solar metallic- 
ity and the stellar mass varying between M = 0.9 Mq and 1.5 Mq, 
with a step of O.IMq. The crosses in each panel indicate log g and 
Tgff error bars of 0.12 dex and 80 K, respectively. 

For each star the obtained effective temperature and 
surface gravity were checked by comparing its position 
in the log g - Teff plane with the theoretical evolution- 
ary track of giv en metallicity and a-enhancement in the 
lYi et alJ ()2004[) grid (Fig.[^. An a-enhancement was 
assumed to be equal to [Mg/Fe] as derived in this study 
and presented in Table[31 Stellar masses corresponding to 
the best fit evolutionary tracks are indicated in Table[Sl 
They range between 0.60 and 0.85 solar mass for the halo 
and thick disk stars that does not contradict with their 
evolutionary status and old age. For the thin disk stars, 
their masses range between 0.85 Mq and 1.5 Mq and ages 
between 1 Gyr and 9 Gyr. 

For two halo cool dwarfs we had to revise / log g, 
to enforce the stars to sit on the main sequence of the ap¬ 
propriate evolutionary track. We stress that the changes 
in atmospheric parameters were not allowed to destroy 
the ionization equilibrium between Fe I and Fe II and the 
Fe I excitation equilibrium. For HD 64090, we adopted 
Teff = 5400 K, which is close to Tirfm = 5440 K (A96), 
and log g = 4.70, which is l.htriogg higher than log gHip- 
For BD -|-66°0268 our final temperature, T^s = 5300 K, 
is close to Tirfm = 5280 K (A96) and log g = 4.72 is 
2 criogg higher than log gmp- 

5.4. Final atmospheric parameters 

The final atmospheric parameters together with the 
NLTE and LTE abundance differences between Fe I and 
Fe II are presented in TableO For the star’s metallicity 
we adopted [Fe/H] determined from lines of Fe H, al¬ 
though the difference in NLTE abundances between Fe I 



[Fe/H] 



Fig. 7.— NLTE (filled circles) and LTE (open circles, top panel 
only) abundance differences Fe I - Fe II for the investigated stars. 
Dashed lines indicate an admissible difference of ±0.06 dex. 
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Fig. 8.— Excitation temperature slopes from Fe I lines for the 
investigated stars. 

and Fe H nowhere exceeds 0.06 dex. The Fe I - Fe H 
abundance differences are also shown in Fig. 0 In NLTE, 
they reveal no trends with any atmospheric parameter, 
either [Fe/H], or Teff, or log g. 

Figure |8| displays the excitation temperature slopes, 
d[Fe/H]/dTexc) calculated from lines of Fe I in the in¬ 
dividual stars. The data reveal no trend with effective 
temperature, and c?[Fe/H]/dTexc = 0.0044 dex/eV, on 
average. The excitation slope seems to depend slightly on 
metallicity at [Fe/H] > — 1, although nowhere its magni¬ 
tude exceeds the error bars. The scatter of data is larger 
for the [Fe/H] < —1.8 than less metal-poor stars due to 
smaller number of observed lines of Fe 1. 

In Fig. 0 we compare surface gravities from the 
Fe I/Fe H NLTE ionization equilibrium and Hipparcos 
parallax methods for the 20 stars, which do not belong 
to the benchmark stellar sample, but have a relative par¬ 
allax uncertainty less than 10 %. The differences are mi¬ 
nor, with Alog g(Sp - Hip) = 0.008 ± 0.037 dex, on av¬ 
erage, and they do not show any trends with surface 
gravity or effective temperature. Having in mind that 
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Fig. 9.— Differences in surface gravity from the Fe I/Fe II NLTE 
ionization equilibrium and Hipparcos parallax methods for 20 stars 
with the spectroscopically derived atmospheric parameters and a 
relative parallax uncertainty less than 10%. In the bottom panel, 
the filled circles and open triangles show stars with log gHip > 4.2 
and <4.2, respectively. 

for each benchmark star the Fe I/Fe II ionization equi¬ 
librium was achieved with log g from the Hipparcos par¬ 
allax method, we infer that the spectroscopic method of 
gravity determination is working in the 5130 K < < 

6600 K and 3.12 < log g < 4.66 range. An exception is 
HD 64090, for which we determined log gsp = 4.70 that 
is 0.13 dex higher than log guip- We, thus, do not con¬ 
firm a temperature trend of Alog g(Sp - Hip) for cool 
(Teff < 5300 K) dw arf (log g > 4.2) stars obtained by 
lllensbv et al.l (l2014ll . 

Final slopes of the [Fe/H] vs. EW plots amount to 
—0.0002 dex/mA, on average. Using the derived micro¬ 
turbulence velocities and basic atmosheric parameters 
Teff, log g, and [Fe/H], we built up the approximation 
formula: 

£,t = -0.21-b0.06x[Fe/H]-F5.6x(Teff/10'‘)-0.43xlog g. 

( 2 ) 

Microturbulence velocity grows towards higher Teff and 
lower surface gravity. The metallicity dependence is 
weak. Similar relatio ns were p reviou s ly found by other 
autho rs , for example iNisse^ (1198^ : lEdvardsson et al. 
( 1993 1: [Allende Prieto et al.l (l2004ll . and iRamfrez et ah 
(|2013f) . In the first three studies theirs formulae did not 
include the [Fe/H] term because these studies dealt with 
the limited metallicity range. Using the formula ([H re¬ 
sults in = l.lkms“^ and 1.4kms“^ for the solar at¬ 
mospheric parameters and 6000/4.0/—1.0, respectively. 
For comparison, the cor responding values are l.lkms”^ 
and 1.5kms“^ from the IRamfrez et aP (l2013ll formula. 

5.5. Notes on individual stars 

We provide additional comments on stars that turned 
out to be interesting for one or the other reason in the 
course of the analysis. 

HD 59984 and HD 105755, with [Fe/H] = -0.69 and 
—0.73, respectively, may be the most metal-poor thin- 
disk stars. Both kinematics, Vpec = 25kms“^ and 
28kms“^, and an age of 8 Gyr for both stars support 
their thin-disk status. 

HD 106516: We failed to obtain self-consistent stellar 
parameters for this star. On the one hand, the star ex¬ 


hibits a typical thick-disk kinematics, Vpec = 108kms“^, 
and chemical signatures, [Mg/Fe] = 0.38 and [Fe/H] = 
—0.73. On the other hand, its age was estimated to be 
6 Gyr and unlikely older than 9 Gyr, identifying the s tar 
as a thin-disk star. According to iGarnev et al.l ( 200111 . it 
is a single-lined spectroscopic binary, with the period P 
= 843.9 days. 

HD 134169: kinematics, Vpec = 5kms“^ identifies it 
as a thin-disk star. However, HD 134169 exhibits a typ¬ 
ical thick-disk chemical signatures, [Mg/Fe] = 0.32 and 
[Fe/H] = —0.78, and old age, r ~ 11 Gyr. Enhance- 
ment of Al re l ative to Fe, with [Al/Fe] = 0.54, found by 
iGehren et al.l (1200611 also suggests a thick-disk origin. 

BD +3T 1458 : this is a halo star. With spectroscopic 
Teff = 5500 K, log g = 3.70, and [Fe/H] = -1.95, the 
star sits on the evolutionary track of M = 0.67Mq at 
r = 24 Gyr. The younger age can only be obtained 
for higher Teff and log g. However, there is no ground 
for a substantial revision of atmospheric parameters. 
The star’s temperat u re is well fixed by sev e ral st udies. 
Alonso et al.l (Il996all: IRamfrez fc MeldndezI (l2005ll . and 


Gonzalez Hernandez fc Bonifaciol ( 20091 1 der ived Ttrfm 
= 551 0 K, 5516 K, and 5582 K, respectively. lAxer et ^ 
(| 199411 determined Teff = 5450 K from the Balmer line 
wing hts. With Teff = 5500 K, log gHip = 3.41 ± 0.18 is 
even lower than the spectroscopically derived value. For 
this star we choose the spectroscopic stellar parameters 
as the final ones. 


5.6. Uncertainties in derived atmospheric parameters 

The following approaches were applied to evaluate the 
uncertainties in the derived atmospheric parameters. 

For the benchmark stars we adopted the Tirfm errors 
as indicated by the original sources. Statistical error of 
the Hipparcos parallax based surface gravity was com¬ 
puted as the quadratic sum of errors of the star’s paral¬ 
lax, effective temperature, mass, visual magnitude, and 
bolometric correction: 

'^logg = (2o’iog7r)^ + (4criogT)^+O’iogMd"(*^-4<^y)^ + (0-4(TBc)^- 

The uncertainties in V magnitude and bolometric correc¬ 
tion together contribute less than 0.01 dex to the total 
log g error. For Teff the contribution nowhere exceeds 
0.03 dex. Stellar masses were well constrained using the 
evolutionary tracks, with an uncertainty of no more than 
0.1 Mq, resulting in 0.04 dex contribution to criogg- 

In Fig.lTUl our final TeffS are compared with the IRFM 
temperatures from Cll and A96. For 47 stars in com¬ 
mon with Gll, our values are 51 K, on average, lower. 

In contrast, this study determined 46 K higher temper¬ 
atures compared with that of A96 for 29 common stars. 
We note very similar statistical errors of 63 K and 69 K 
for the temperature differences (this work - Cll) and 
(this work - A96), respectively. Based on these compar¬ 
isons, we estimate the systematic and statistical errors 
of Teff for the 31 stars with the spectroscopically derived 
atmospheric parameters to be 50 K and 70 K. 

For 20 stars with the spectroscopic surface gravities, 
log gsp , their Hipparcos parallaxes were measured with 
a relative parallax error less than 10 %, and we calculated 
reliable log gHip values. As shown in Fig. [HI log gsp - 
log gHip = 0.008 ±0.037 dex. This led us to estimate the 
uncertainty in log gsp to be 0.04 dex. 
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Fig. 10. — Differences in effective temperature between this 
study and the literature dat a. The top pane l correspond s to inter- 
ferometric measurements o f IBovaiian et al. , sq uares), 

ICreevev et al.l ll2Q12l. 120151. filled circles). iNorth et al.l l|2Q09 . trian- 
gle). and Ivon BrauneT^n (120141 . cross). See sect. 0 for more de¬ 
tails. The_jnldffle_and_bottom pane ls use the IRFM temp eratures 
from ICasagrande et ^ ] (|2011l ') and lAlonso et al.l l|1996bli , respec¬ 
tively. 


Statistical error of [Fe/H] was defined by the disper¬ 
sion, cr, for lines of Fe II in given star. 

Statistical error of the microturbulence velocity was 
adopted to be common for the whole stellar sample, and 
it was defined by the dispersion in the single measure¬ 
ments about the relation ([2]). It amounts to 0.14 kms“^. 

It is worth noting that the uncertainty in spectro¬ 
scopically derived atmospheric parameters can be larger 
than that quoted above. For example, for the halo star 
HD 74000 we found the two sets of stellar parameters 
that reproduce well the Fe I/Fe II ionization and Fe I 
excitaion equilibrium iFig. lllIl and fit the appropriate 
evolutionary tracks. These are Teff = 6225 K, log g = 
4.13, and [Fe/H] = —1.97 suggesting the star’s age of 
15 Gyr and Teg = 6360 K, log g = 4.30, and [Fe/H] = 
— 1.93 resulting in r = 13.5 Gyr. The first set relied, in 
fact, on riRFM(A96) and the second one on Tirfm(GII). 
A shift of -1-135 K in Tgg leads to -1-0.17 dex shift in 
log g, and both sets of parameters provide an acceptable 
star’s age. No of two sets can be preferred from a com¬ 
parison of log gsp with the corresponding Hipparcos 
parallax based surface gravity, log gHip = 4.16 ±0.16 
for Teff = 6225 K and log gHip = 4.20 ± 0.16 for Teg 
= 6360 K. From analys is of Hq, and H^ in HD 74000 
iMashonkina et ahl (1200311 inferred T^s = 6225 K. Two 
pairs of Teg / log g, which provide the Fe I/Fe H ion¬ 
ization equilibrium and reasonable star’s mass and age, 
were also obtained for HD 108177 (6100/4.22/—1.67 and 
6250/4.40/-1.62), and G090-003 (6100/3.90/-2.04 and 
5930/3.80/—2.10). We, therefore, recommend to ap¬ 
ply various spectroscopic and non-spectroscopic meth¬ 
ods to given star to find an unique solution for the star’s 
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Fig. 11.— Impact of changes in Teff and log g on the NLTE 
differential abundances derived from the Fe I (open circles) and 
Fe II (filled circles) lines in HD 74000. The top panel corresponds 
to our final parameters, Tgff = 6225 K, log g = 4.13, [Fe/H] = 
— 1.97, and the bottom one to Teff = 6360 K, log g = 4.30, [Fe/H] 
= —1.93. The differences in NLTE abundances Fe I - Fe II = 
—0.03 dex and H-O.Ol dex, respectively. 

Teff / log g. 

6. COMPARISON WITH OTHER STUDIES 

For ten stars of our sample, their Tegs were determined 
in the literature based on measurements of the angu¬ 
lar dia meters, trigonometric p arallaxes, and bolometric 
fluxes. IBovaiian et al.l (120131 B13) published the most 
numerous list of the interferometric temperatures. Tint, 
based on angular diameters measured with the CHARA 
array (330 m maximum baseline). FigurejTU] displays 
the temperature differences betwee n this study and in ¬ 
terferometric measuremen ts of B13.lNor th et al.) (2QQ^; 
ICreevev et al.l (|2012L[20l^ . and Ivon Braun et"^ ( 2014li . 
Our values are, on average, higher, with ATeg = 135 ± 
126 K for ten stars and 78 ± 81 K, when excluding the 
two outliers, HD 103095 and HD 140283. We selected 
three stars for a detailed comparison. 

HD 102870. _ Two successive determinations of 

IBovaiian et al.l (|2012l l and B13 resulted in Tint = 6132 K 
and 6054 K, with small temperature errors of 36 K and 
13 K, respectively. A downward revision by 78 K was due 
to employing the different bolometric fluxes. The most 
re cent temperature o f B13 is in line with the earlier data 
of iNorth et al.l (|2009f ). Tint = 6059 ± 49 K. Surface grav¬ 
ity based on the asterose ismic measurements amounts to 
log gseis = 4.11 ± 0.02 (ICreevev et al.l[201^ . With Teg 
= 6060 K and log g = 4.11, we failed to achieve the 
ionization equilibrium between Fe I and Fe H and ob¬ 
tained Fe I - Fe H = —0.12 dex in NLTE. It is worth 
noting, log gHip = 4.14 adopted in this study agrees well 
with log gseis- For solar-type stars the asteroseismology 
method depends only weakly on Tgg. Indeed, a 200 K 
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Fig. 12.— Theoretical NLTE flux profiles of Ha calculated with 
^eff = 4820 K (dashed curve) and T’eff = 5130 K (continuous curve) 
compared to the observed FOCES spectrum of HD 103095 (bold 
dots). The differences between observed and calculated spectra, 
(O - C), are shown in the upper part of the panel. 


Fig. 13.— Theoretical NLTE flux profiles of Ha calculated with 
^eff = 5535 K (dash-dotted curve), 5645 K (dashed curve), and 
5780 K (continuous curve) compared to the observed VLT/UVES 
spectrum of HD 140283 (bold dots). The (O - C) values are shown 
in the upper part of the panel. 


temperature difference produces a difference of 0.007 dex 
in log g. To obtain consistent NLTE abundances from 
lines of Fe I and Fe II using log g = 4.14, one needs to 
have 110 K higher temperature of HD 102870 (TableE]) 
compared with its Tint ■ 

HD 10 3095. Using th e CHA RA array measurements, 
B13 and iCreevev et al.l (| 20 1211 derived Tint = 4771 ± 
18 K and 4818 ± 54 K, respectively. We checked Teg = 
4820 K and log gHip = 4.60 with various spectroscopic 
temperature and gravity indicators. 

(i) The Fe I/Fe II ionization equilibrium is not fulfilled, 
and the NLTE abundance difference amounts to Fe I - 
Fe II = -0.27 dex. 

(ii) We compared the NLTE abundances from C I 9094 
A, 9111 A with the carbon abundance from a number of 
molecular CH bands, which are known to be sensitive to 
Tes variation. The NLTE meth od for C I and atomic and 
molec ular data were taken from [Alexeeva &: Mashonkinal 
(|2015ll . With the 4820/4.60/—1.3 model, the abundance 
difference. Cl- CH = 0.56 dex, is very large, while it 
amounts to —0.06 dex for our final parameters, Teff = 
5130 K and log g = 4.66. 

(iii) Effective temperature can also be constrained 
from analysis of the H^ line wings iFig. lT^ . The the- 
oretical NLTE profiles o f Hq, were computed following 
iMashonkina et ahl (|2008ll . It can be seen that Teg = 
4820 K leads to shallower wings of H^ compared with 
the observations, but Teg = 5130 K to deeper core-to- 
wing transition region. The bes t fit was achieved for Teg 
= 5030 K. It is worth noting, ICavrel et al.l (|2011ll ob¬ 
tained 100 K, on average, lower temperatures from the 
H„ line wings compared with accurate effective tempera¬ 
tures from the apparent angular diameter for the eleven 
FGK-type stars in the —0.7 < [Fe/H] < 0.2 range. In 
contrast. Tint = 4771 K and 4818 K of HD 103095, as de¬ 
termined by B13 and ICreevev et all (I2012D . respectively, 
are more than 200 K lower than the temperature from 
H«. 

Thus, no spectroscopic temperature indicator supports 
the literature data on Tint fo r HD 1 03095. 

HD 140283. ICreevev et all (l2nT^ derived T„t = 5534 
± 103 K and 5647 ± 105 K assuming zero-reddening and 
Ay = 0.1™, respectively. It is worth noting, HD 140283 
is a nearby star, at a distance of 58 pc from the Sun, and 
an interstellar absorption of Ay = 0.1™ is unlikely pro¬ 
duced. Of the two temperatures Tint = 5534 K should be 
preferred, nevertheless we checked both using the same 


spectroscopic indicators as for HD 103095. The NLTE 
abundance differences Fe I - Fe H were found to be 
—0.18 dex and —0.09 dex, respectively. For both temper¬ 
atures an abundance difference between carbon atomic 
and molecular lines is large, with (Cl - CH) = 0.69 dex 
and 0.32 dex, while consistent abundances, (Cl - CH) 
= 0.00 dex, were obtained for our final parameters. Teg 
= 5780 K and log g = 3.70. FigurelT^ shows theoretical 
profiles of Hq, in the three different model atmospheres 
compared with the observed spectrum of HD 140283. It 
is evident that Tint = 5534 K is too low and it does not 
fit any spectroscopic indicator of Teff- 

It is worth noting, both HD 103095 and HD 140283 
have rathe r small anguly dia meters, dmt = 0.679 ± 
0.007 mas (ICreevev et al.l[ 2012 ll and 0.353 ± 0.013 mas 
(|Creevev et al.ll2015f) . respectively, which can be over- 
estimat ed resulting in u n deres timated temperatures. 
Indeed, ICasagrande et al.l (l20I4ll su spected systematic 
trends in the Bovaiian et al.l ( ^121 BI2) dataset that 
was also based on the CHARA array interferometric mea¬ 
surements. They obtained the differences (Tirfm(CII) 
- Tijit(B 12 )) growing towards smaller angular diameter 
and reaching -1-200 K, on average, at Oint = 0.8 mas. 

Asteroseismic measurements of surface gravity were 
made for one more star of our sample, HD 49933 
(|Kallinger et al.|[2010ll . They resulted in log gseis = 4.22, 
which is only 1.4(Tg higher, than our value. 

Spectroscopic stellar parameters were determined in 
a number of rec ent studies. We sele cted three of them 
for comparison. iBensbv et al.l (|2014il used lines of Fe I 
and Fe H to derive temperatures and gravities of the 
extended stellar sample. The LTE abundances from 
individual lines were corrected usi ng the grid of th e 
NLTE abundance corrections from iLind et ^ (j2012t l. 
For nearby stars wit h very good Hipparcos parallaxes, 
iBensbv et al.l (|20I4ll calculated also the log gntp values 
and found systematic discrepancies between the surface 
gravities from the two methods for the stars with log g 
> 4.2 and Teg < 5650 K. Therefore, they applied em¬ 
pirical corrections to the atmospheric parameters from 
ionization balance. It should be reminded, we obtained 
no temperature trend of Alog g(Sp - Hip) for our pro¬ 
gram stars (Sect.E]). FigurefHl (top panel) shows the 
differences in Teg an d log g for 17 stars in common with 
IBensbv et al.l (j2014l l. On average, ATeff = —34±87 K, 
Alog g = -0.02±0.11, and A[Fe/H] = -0.02±0.11 in the 
sense “this work minus Bensby”. The two independent 
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Fig. 14. — Differences in effective temperature and surface 
gr avity between this st udy and the spectroscopic determ inations 
bv IBens by et all 11201-^) (top panel ) , ITakeda et al.l (120051) (middle 
panell. and lFuhrmannf (llOok . f2004l'l (bottom panel). 

research provide consistent stellar parameters, with the 
small systematic shifts and the statistical errors being 
typical of stellar parameter determinations. However, 
for some individual stars the discrepancies are uncom¬ 
fortably large. For example, ATeff exceeds 100 K for 
five stars and Alog g > 0.2 for HD 22879. We note, in 
particular, the differences of different sign for two well- 
studied very metal-poor stars HD 84937 and HD 140283, 
with ATeff = —190 K and -1-120 K and Alog g = —0.14 
and -1-0.13, respectively. This can be due to large un¬ 
certainty in evaluating a slope of the log A(Fe I) versus 
Texc trend, where the total number of investigated lines 
is limited. Since the Fe I/Fe H ionization equilibruim de¬ 
pends on not only surface gravity, but also Teff, an over¬ 
estimation (underestimation) of Teg produces an upward 
( downward) shi f t in lo g g. 

iTakeda et al.l (|2005ll derived stellar Teg and log g from 
lines of Fe I and Fe H under the LTE assumption. How¬ 
ever, all the 18 stars in common with our work lie in the 
— 1.3 < [Fe/H] < 0.2 metallicity range, where the depar¬ 
tures from LTE for Fe I and Fe H are small. Therefore, 


the obtained differences in Tgg and log g (Fig. [Ml mid¬ 
dle panel) cannot be caused by applying different line- 
formation treatments. The data in Fig. [TTl suggest, most 
probably, that the uncertainty in the derived effective 
temperature was directly translated to the uncertainty 
i n surface gravity. 

iFuhrmaiml (I1998L l2004[) applied different spectroscopic 
approaches to derive Tgg from the Balmer line wing fits 
and log g from the Mg Ib line profile fits under the LTE 
assumption. For 16 stars in common with this work, 
the mean differences in effective temperature and iron 
abundance, ATgg = 29 ± 71 K (Fig. [Til bottom panel) 
and A[Fe/H] = 0.00 ± 0.07, do not exceed the error bars. 
The only star, HD 30743, was found with ATgg > 100 K. 
However, this study obtains higher surface gravities, with 
Alog g = 0.08±0.07. The discrepancy can be exp l ained , 
inmrt, by using the LTE approach in iFuhrmannI (Il998l 
l2004f) . In the stellar parameter range, with which we 
concern. Mg I is subject to overionization resulting in 
weakened Mg Ib line. Ignoring the departures from LTE 
leads to an underestimation of the derived gravity. 

7. CONCLUSIONS 

The sample of 51 nearby FG dwarf stars uniformly 
distributed over the —2.60 < [Fe/H] < -1-0.20 metallicity 
range was selected for a systematic NLTE investigation 
of the Galaxy chemical evolution. A membership of in¬ 
dividual stars to the particular galactic stellar popula¬ 
tion, namely the thin disk, the thick disk, and the halo, 
was mainly identified using the star’s kinematics and for 
few stars using also the chemical signatures, [Fe/H] and 
[Mg/Fe]. A combination of the photometric and spectro¬ 
scopic methods was applied to derive a homogeneous set 
of the stellar atmosphere parameters: Tgg, log g, [Fe/H], 
and ■ Our spectroscopic analyses took advantage of em¬ 
ploying the high-resolution {R > 60 000) observed spec¬ 
tra and NLTE line formation for Fe I and Fe H in the 
classical ID model atmospheres. Spectroscopic method 
of Tgg/log g determination was tested using the 20 bench¬ 
mark stars, for which there are multiple measurements 
of the IRFM effective temperature and their Hipparcos 
parallax error is less than 10%. An efficiency of poorly 
known inelastic Fe-|-H collisions in the SE of Fe I-Fe H 
was estimated empirically from analysis of their different 
influence on lines of Fe I and Fe H in the MP bench¬ 
mark stars. We found abundances from the two ioni¬ 
sation stages to be consistent within 0.06 dex for every 
star, when applying S'h = 0.5. 

Stellar parameters of the remaining 31 stars were ob¬ 
tained spectroscopically from the NLTE analysis of the 
iron lines. For lines of Fe H in our program stars the 
NLTE abundance corrections do not exceed 0.01 dex in 
absolute value. The deviations from LTE for Fe I grow 
towards higher Tgg and lower [Fe/H] and log g. For stars 
with either [Fe/H] > —0.75, or Tgg < 5750 K, or log g > 
4.20, the difference in average abundance between NLTE 
and LTE was found to be less than 0.06 dex, which trans¬ 
lates to a shift of less than 0.1 dex in log g. Since NLTE 
leads to weakened lines of Fe I, but it does not affect lines 
of Fe H until the extremely low metallicities, the ioniza¬ 
tion equilibrium Fe I/Fe H is achieved for higher gravities 
in NLTE than in LTE. The NLTE analysis is crucial for 
the VMP turn-off and subgiant stars. Indeed, the shift 
in log g between NLTE and LTE reaches -1-0.45 dex for 
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BD-13° 3442, with Tes = 6400 K, log g(NLTE) = 3.95, 
and [Fe/H] = -2.62. 

The obtained effective temperatures and surface gravi¬ 
ties were checked by comparing the star’s position in the 
log g - Teff plane with the theoretical evolution ary track 
of giy en metallicity and a-enhancement in the lYi et all 
(|2004ll grid. Most stars reveal self-consistent data on 
star’s mass, age, and kinematics. 

Our final effective temperatures lie exactly in between 
the Tirfm scales of A96 and Cll, with a systematic shift 
of -1-46 K and —51 K, respectively. We estimate the 
Teff statistical error to be about 70 K. Surface gravi¬ 
ties obtained from the Fe I/Fe II and Hipparcos parallax 
methods were found to agree well. We do not confirm 
a temperature tren d of Alog g(Sp - Hip) reported by 
iBensbv et al.l (j20I4l l for the log g > 4.2 and < 5650 K 
stars. 

We recommend to apply a line-by-line differential ap¬ 
proach relative to the Sun, to minimise the effect of the 
uncertainty in ^/-values on the final results. It is quite 
efficient for the [Fe/H]> —1.5 stars. However, in the 
more metal-poor stars, the line-to-line scatter for Fe I 
has a similar magnitude for the differential and absolute 
abundances. This is, probably, due to uncertainties in 
Ce-values for the high-excitation lines. We note, in par¬ 
ticular, Fe I 5367A (Fexc = 4.41 eV) and Fe I 5383 A 
(Fexc = 4.31 eV). 

It was found that none of checked spectroscopic tem¬ 
perature indicators, such as Fe I versus Fe H, C I versus 
CH, and the Hq, line wings, do support interferometric 
effecti ve temperatures of HD 103095, Tint = 4818 K ± 
54 K (ICreevev et ahllml. an d HD 140283, T,nt = 5534 
± 103 K ( Creevev et al.ll2015f) . Both stars have rather 
small angular diameters, and their measurements can 
be affected by some syst ematic t rends, as suspected by 
ICasagrande et al.l (|2014l l for the iBovaiian et al.l (1201^ 
dataset. 

We conclude that the NLTE analysis of lines of iron 
in the two ionisation stages, Fe I and Fe H, is efficient 
in deriving reliable atmospheric parameters for the FG- 
type dwarf stars in the broad metallicity range, down 
to [Fe/H] = —2.6. The stellar parameter determina¬ 
tions would benefit from using the complementary data 
on interferometric and IRFM temperatures, trigonomet¬ 
ric parallaxes, and asteroseismology measurements. We 
also recommend to combine different spectroscopic tem¬ 
perature and gravity indicators and check the obtained 
Teff/log g values with the theoretical evolutionary tracks 
/ isochrones. 

The obtained accurate atmospheric parameters will be 
used in the forthcoming papers to determine NLTE abun¬ 
dances of important astrophysical elements from lithium 
to europium and improve observational constraints on 
the chemo-dynamical models of the Galaxy evolution. 
We also plan to extend our sample toward lower metal¬ 
licity and surface gravity. 
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TABLE 2 

Characteristics of observed spectra. 


HD,BD 

V 

[Fe/H] 

S/N 

texpi S 

Nspec 

Date (yyyy-mm-dd) 




Shane/Hamilton 




19373 

4.05 

0.10 

234, 228 

30.00, 30.00 

2 

2012-01-09, 2012-01-09 

22484 

4.28 

0.01 

171, 160 

30.00, 30.00 

2 

2012-01-09, 2012-01-09 

22879 

6.74 

-0.84 

170, 195 

360.00, 360.00 

2 

2012-01-09, 2012-01-09 

24289 

9.96 

-1.94 

112, 128 

2700.00, 2700.00 

2 

2012-01-09, 2012-01-09 

30562 

5.77 

0.17 

180, 190 

120.00, 120.00 

2 

2012-01-09, 2012-01-09 

30743 

6.26 

-0.44 

188, 184 

240.00, 240.00 

2 

2012-01-09, 2012-01-09 

34411 

4.70 

0.01 

220, 218, 209+1 

45.00, 45.00, 84.19 

3 

2012-01-11, 2012-01-11, 2012-01-07 

43318 

5.65 

-0.19 

210, 212, 183F 

100.00, 100.00, 215.14 

3 

2012-01-11, 2012-01-11, 2012-01-07 

45067 

5.90 

-0.16 

230, 235 

180.00, 180.00 

2 

2012-01-09, 2012-01-09 

45205 

8.00 

-0.87 

140, 177 

1200.00, 1800.00 

2 

2012-01-09, 2012-01-09 

49933 

5.78 

-0.47 

216, 203 

659.71, 135.0 

2 

2012-01-08, 2012-01-10 

52711 

5.93 

-0.21 

256, 163F 

160.00, 119.18 

2 

2012-01-10,2012-01-06 

58855 

5.36 

-0.29 

224, 167F 

150.00, 75.13 

2 

2012-01-09, 2012-01-06 

59374 

8.50 

-0.88 

122F, 131F 

1174.47, 812.31 

2 

2012-01-07, 2012-01-06 

59984 

5.93 

-0.69 

182F, 165F 

269.17, 229.17 

2 

2012-01-07, 2012-01-06 

62301 

6.75 

-0.70 

197, 158F 

450.00, 238.20 

2 

2012-01-09, 2012-01-06 

64090 

8.30 

-1.73 

280, 136F 

2500.00, 687.35 

2 

2012-12-26, 2012-01-06 

69897 

5.10 

-0.25 

186, 180F 

240.00, 66.18 

2 

2012-01-08, 2012-01-06 

74000 

9.67 

-1.97 

144, 142, 74F 

2400.00, 2400.00, 1454.51 

3 

2012-01-10, 2012-01-10, 2012-01-06 

76932 

5.86 

-0.98 

170, 181F, 210F 

612.32, 258.23, 489.24 

3 

2012-01-08, 2012-01-06, 2012-01-07 

82943 

6.54 

0.19 

214, 170F 

270.00, 342.25 

2 

2012-01-10, 2012-01-06 

84937 

8.28 

-2.16 

122, 153, 95 

1475.81, 1800.0, 846.41 

3 

2012-01-08, 2012-01-09, 2012-01-06 

89744 

5.74 

0.13 

220, 165F 

135.00, 121.16 

2 

2012-01-10, 2012-01-06 

90839 

4.83 

-0.18 

190, 169F 

60.00, 137.20 

2 

2012-01-10, 2012-01-06 

92855 

7.28 

-0.12 

171, 163F 

600.00, 309.20 

2 

2012-01-10, 2012-01-06 

94028 

8.23 

-1.47 

118, 172 

1800.00, 1317.69 

2 

2011-03-15, 2011-03-15 

99984 

5.95 

-0.38 

181F, 183F 

517.22, 491.24 

2 

2012-01-07, 2012-01-07 

100563 

5.70 

0.06 

171F, 173F 

413.25, 383.21 

2 

2012-01-07, 2012-01-07 

102870 

3.61 

0.11 

208F, 211F 

77.12, 91.13 

2 

2012-01-07, 2012-01-07 

103095 

6.45 

-1.26 

188 

300.00 

1 

2012-01-10 

105755 

8.59 

-0.73 

175, 103F 

2400.0, 1052.44 

2 

2012-01-10, 2012-01-06 

106516 

6.11 

-0.73 

151F, 155F 

434.25, 546.25 

2 

2012-01-07, 2012-01-07 

108177 

9.66 

-1.67 

60, 111 

3000.00, 3600.00 

2 

2011-03-15, 2011-03-15 

110897 

6.00 

-0.57 

260, 170F, 172F 

600.00, 568.30, 508.23 

3 

2012-01-08, 2012-01-07, 2012-01-07 

114710 

4.26 

0.06 

184F, 186F 

112.13, 129.15 

2 

2012-01-07, 2012-01-07 

115617 

4.74 

-0.10 

177, 170 

160.00, 160.00 

2 

2012-01-08, 2012-01-08 

134088 

8.00 

-0.80 

130 

1800.00 

1 

2011-03-15 

134169 

7.67 

-0.78 

214, 212 

750.00, 750.00 

2 

2012-01-10, 2012-01-10 

138776 

8.72 

0.24 

211 

1800.00 

1 

2012-01-11 

142091 

4.82 

-0.07 

277, 275 

60.00, 60.00 

2 

2012-01-11, 2012-01-11 

142373 

4.62 

-0.54 

140, 147 

240.00, 240.00 

2 

2012-01-08, 2012-01-08 

-04°3208 

9.99 

-2.20 

60, 66 

1500.00, 2400.00 

2 

2012-01-08, 2012-01-08 

-13°3442 

10.37 

-2.62 

103, 102, 104, 102 

2700.00, 2700.00, 

4 

2012-01-09, 2012-01-09, 





2700.00, 2700.00 


2012-01-09, 2012-01-09 

+24° 1676 

10.80 

-2.44 

86, 90, 92, 92, 57F 

2400.00, 2400.00, 2400.00, 

5 

2012-01-11, 2012-01-11, 2012-01-11, 





2400.00, 1581.60 


2012-01-11, 2012-01-06 

+29°2091 

10.22 

-1.91 

82, 121, 115, 120, 83F 

2400.00, 2700.00, 2700.00, 

5 

2012-01-08, 2012-01-10, 2012-01-10, 





2700.00, 1287.46 


2012-01-10, 2012-01-06 

+37° 1458 

8.92 

-1.95 

235, 105F 

2500.00, 897.38 

2 

2012-01-10, 2012-01-06 

+66°0268 

9.88 

-2.06 

112, 110, 115 

2400.00, 2400.00, 2400.00 

3 

2012-01-10, 2012-01-10, 2012-01-10 

G090-003 

10.50 

-2.04 

106, 120, 122, 66F 

2700.00, 2700.00, 

4 

2012-01-11, 2012-01-11, 





2700.00, 1181.44 


2012-01-11, 2012-01-06 




CFHT/ESPaDOnS 




+07°4841 

10.38 

-1.46 

152 

1470 

2 


+09°0352 

10.17 

-2.09 

160 

FOCES 

2400 

2 


22879 

6.74 

-0.84 

200 

900.00, 1800.00, 1800.00 

3 

1996-09-02, 1999-01-02, 1996-10-31 

34411 

4.70 

0.01 

200 

1500.00, 1500.00 

2 

1999-08-21, 1998-12-27 

84937 

8.28 

-2.16 

200 

3600.00, 3300.00 

2 

1999-03-01, 1999-03-01 

94028 

8.23 

-1.47 

200 

2400.00, 2400.00 

2 

2000-01-17, 2000-01-17 

103095 

6.45 

-1.26 

200 

900.00 

1 

2000-05-19 

142373 

4.62 

-0.54 

200 

Others 

600.00, 600.00 

2 

1998-06-09, 1998-06-09 

49933^ 

4.15 

-0.47 

500 




140283^ 

7.21 

-2.46 

200 




-04°32083 

9.99 

-2.20 

200 





Note. 


F indicates the observations carried out by Debra Fischer, ^ 3.6-m/HARPS, ^ VLT/UVES. 
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TABLE 3 

Stellar kinematics and membership to particular galactic stellar population. 


HD, BD 

X-8000 

(pc) 

Y 

(pc) 

Z 

(pc) 

u 

(km s~^) 

V 

(kms~^) 

w 

(kms~^) 

[Mg/Fe] 

NLTE 

Stellar 

population 

19373 

10 

-10 

-1 

-75.2 

-16.1 

21.5 

0.00 

± 

0.03 

Thin disk 

22484 

10 

0 

-9 

1.7 

-15.5 

-42.3 

-0.05 

± 

0.02 

Thin disk 

22879 

20 

0 

-17 

-109.9 

-88.2 

-41.3 

0.30 

± 

0.02 

Thick disk 

24289 

160 

30 

-139 

-116.7 

-172.6 

173.9 

0.21 

± 

0.04 

Halo 

30562 

20 

10 

-13 

-54.9 

-72.0 

-21.6 

0.04 

± 

0.05 

Thin disk 

30743 

20 

20 

-18 

24.8 

-5.2 

-22.5 

0.08 

± 

0.06 

Thin disk 

34411 

10 

0 

0 

-74.5 

-35.4 

4.4 

-0.03 

± 

0.03 

Thin disk 

43318 

30 

20 

-5 

47.0 

1.2 

-37.0 

-0.05 

± 

0.05 

Thin disk 

45067 

30 

20 

-3 

-16.5 

-64.4 

13.7 

0.01 

± 

0.04 

Thin disk 

45205 

60 

-30 

26 

-35.9 

-89.5 

17.0 

0.39 

± 

0.03 

Thick disk 

49933 

30 

20 

0 

25.7 

-13.1 

-9.2 

0.16 

± 

0.10 

Thin disk 

52711 

20 

0 

5 

-18.6 

-77.8 

-9.0 

0.01 

± 

0.05 

Thin disk 

58855 

20 

0 

9 

25.3 

-15.4 

-3.9 

0.07 

± 

0.03 

Thin disk 

59374 

40 

20 

14 

-48.5 

-120.3 

-9.8 

0.30 

± 

0.02 

Thick disk 

59984 

20 

20 

3 

8.0 

-11.4 

-21.2 

0.10 

± 

0.02 

Thin disk 

62301 

30 

0 

16 

-7.4 

-110.0 

-22.1 

0.26 

± 

0.03 

Thick disk 

64090 

30 

0 

13 

266.9 

-227.6 

-89.6 

0.24 

± 

0.04 

Halo 

69897 

20 

0 

9 

-23.9 

-38.7 

7.1 

-0.01 

± 

0.06 

Thin disk 

74000 

60 

110 

34 

216.8 

-348.0 

62.2 

0.37 

± 

0.03 

Halo 

76932 

10 

20 

7 

-49.1 

-91.3 

70.0 

0.32 

± 

0.05 

Thick disk 

82943 

10 

20 

13 

13.3 

-13.5 

-12.6 

-0.05 

± 

0.05 

Thin disk 

84937 

40 

30 

52 

205.0 

-214.7 

-7.9 

0.21 

± 

0.08 

Halo 

89744 

20 

0 

33 

-10.6 

-29.9 

-14.3 

-0.02 

± 

0.04 

Thin disk 

90839 

10 

0 

10 

-13.9 

-1.8 

2.4 

0.02 

± 

0.03 

Thin disk 

92855 

20 

0 

32 

-42.1 

-22.4 

-13.9 

-0.02 

± 

0.04 

Thin disk 

94028 

20 

10 

42 

-33.7 

-129.3 

13.5 

0.35 

± 

0.05 

Thick disk 

99984 

20 

-10 

49 

-6.4 

11.0 

-36.1 

0.02 

± 

0.03 

Thin disk 

100563 

0 

10 

24 

-14.5 

-20.9 

-9.8 

-0.03 

± 

0.03 

Thin disk 

102870 

0 

10 

10 

40.4 

3.2 

7.3 

-0.02 

± 

0.03 

Thin disk 

103095 

0 

0 

9 

278.3 

-157.2 

-14.8 

0.24 

± 

0.06 

Halo 

105755 

30 

-30 

78 

26.7 

5.3 

8.3 

0.34 

± 

0.08 

Thin disk 

106516 

0 

10 

18 

53.4 

-73.4 

-57.7 

0.38 

± 

0.04 

Thick disk 

108177 

-20 

50 

99 

126.7 

-262.9 

32.0 

0.23 

± 

0.04 

Halo 

110897 

0 

0 

17 

-41.7 

7.0 

75.3 

0.18 

± 

0.03 

Thin disk 

114710 

0 

0 

9 

-49.9 

11.5 

8.4 

-0.07 

± 

0.04 

Thin disk 

115617 

0 

0 

6 

-23.3 

-47.8 

-31.3 

-0.10 

± 

0.02 

Thin disk 

134088 

-30 

0 

26 

-25.6 

-72.4 

-68.7 

0.36 

± 

0.03 

Thick disk 

134169 

-40 

0 

44 

2.5 

-3.3 

-1.9 

0.32 

± 

0.02 

Thin disk 

138776 

-50 

0 

45 

6.4 

-56.4 

-5.0 

0.02 

± 

0.04 

Thin disk 

140283 

-50 

0 

32 

-249.0 

-253.3 

41.1 

0.23 

± 

0.03 

Halo 

142091 

-10 

-20 

24 

33.6 

-40.7 

-18.1 

0.07 

± 

0.05 

Thin disk 

142373 

0 

-10 

12 

-41.6 

10.8 

-68.3 

0.18 

± 

0.01 

Thin disk 

-04° 3208 

-30 

110 

169 

-123.5 

-309.9 

-120.2 

0.21 

± 

0.06 

Halo 

-13° 3442 







0.23 

± 

0.01 

Halo 

+07° 4841 

-50 

-150 

-123 

-224.4 

-283.9 

-55.3 

0.43 

± 

0.05 

Halo 

+09° 0352 

120 

-40 

-120 

-156.5 

-176.9 

97.7 

0.45 

± 

0.05 

Halo 

+24° 1676 

360 

90 

126 

166.8 

-483.7 

106.4 

0.21 

± 

0.06 

Halo 

+29° 2091 

40 

20 

85 

157.3 

-345.9 

88.4 

0.34 

± 

0.04 

Halo 

+37° 1458 

150 

-10 

26 

-280.3 

-233.0 

-26.0 

0.38 

± 

0.05 

Halo 

+66° 0268 

40 

-40 

10 

-181.0 

-426.4 

-73.1 

0.15 

± 

0.07 

Halo 

G090-003 

820 

90 

332 

10.4 

-824.6 

-192.4 

0.25 

± 

0.05 

Halo 
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TABLE 4 

Line data, iron NLTE and LTE abundances from an 

ANALYSIS OF THE SOLAR SPECTRUM, AND EQUVALENT 
WIDTHS {EW) OF THE SOLAR LINES. 


A, A 

F, 

-^exc 

(eV) 

log gf 

log Ce 

log A 
NLTE 

log A 
LTE 

EW 

mA 

Fe I 

4920.50 

2.83 

0.07 

-30.51 

-4.61 

-4.61 

466.0 

5198.72 

2.22 

-2.14 

-31.32 

-4.52 

-4.53 

102.0 

5217.40 

3.21 

-1.07 

-30.37 

-4.57 

-4.57 

130.3 

5232.94 

2.94 

-0.06 

-30.54 

-4.68 

-4.68 

371.1 

5236.20 

4.19 

-1.50 

-31.32 

-4.68 

-4.69 

33.7 

5242.50 

3.63 

-0.97 

-31.56 

-4.47 

-4.47 

93.4 

5281.79 

3.04 

-0.83 

-30.53 

-4.68 

-4.69 

163.6 

5324.18 

3.21 

-0.10 

-30.42 

-4.54 

-4.54 

332.5 

5367.47 

4.41 

0.44 

-30.20 

-4.78 

-4.78 

170.1 

5379.58 

3.69 

-1.51 

-31.56 

-4.46 

-4.47 

64.1 

5383.37 

4.31 

0.64 

-30.37 

-4.72 

-4.72 

220.5 

5393.17 

3.24 

-0.72 

-30.42 

-4.60 

-4.60 

171.2 

5491.83 

4.19 

-2.19 

-31.33 

-4.53 

-4.54 

14.5 

5586.76 

3.37 

-0.10 

-30.38 

-4.58 

-4.58 

289.3 

5662.52 

4.18 

-0.57 

-30.52 

-4.41 

-4.42 

105.7 

5696.09 

4.55 

-1.72 

-30.21 

-4.68 

-4.69 

14.6 

5705.46 

4.30 

-1.36 

-30.47 

-4.61 

-4.62 

41.2 

5778.45 

2.59 

-3.44 

-31.37 

-4.60 

-4.61 

22.7 

5855.08 

4.61 

-1.48 

-30.21 

-4.58 

-4.59 

24.4 

5916.25 

2.45 

-2.99 

-31.45 

-4.43 

-4.44 

56.7 

6065.49 

2.61 

-1.53 

-31.41 

-4.51 

-4.52 

132.1 

6082.71 

2.22 

-3.57 

-31.74 

-4.55 

-4.56 

34.7 

6151.62 

2.18 

-3.30 

-31.58 

-4.52 

-4.53 

51.1 

6200.32 

2.61 

-2.44 

-31.43 

-4.47 

-4.48 

74.6 

6213.43 

2.22 

-2.48 

-31.58 

-4.58 

-4.59 

86.0 

6229.23 

2.85 

-2.80 

-31.32 

-4.65 

-4.66 

37.8 

6252.55 

2.40 

-1.69 

-31.52 

-4.52 

-4.53 

135.4 

6393.61 

2.43 

-1.43 

-31.53 

-4.63 

-4.64 

150.8 

6411.65 

3.65 

-0.60 

-30.38 

-4.57 

-4.57 

154.7 

6421.35 

2.28 

-2.03 

-31.80 

-4.56 

-4.57 

112.7 

6518.37 

2.83 

-2.46 

-31.37 

-4.59 

-4.60 

64.0 

Fe II 

4233.17 

2.58 

-1.97 

-32.01 

-4.64 

-4.64 

132.1 

4508.29 

2.84 

-2.44 

-32.00 

-4.48 

-4.48 

93.0 

4582.83 

2.83 

-3.18 

-32.03 

-4.54 

-4.54 

57.7 

4620.52 

2.82 

-3.21 

-32.02 

-4.63 

-4.63 

53.9 

4923.93 

2.88 

-1.26 

-32.03 

-4.69 

-4.69 

196.1 

5018.44 

2.88 

-1.10 

-32.04 

-4.72 

-4.72 

219.1 

5197.58 

3.22 

-2.22 

-32.02 

-4.58 

-4.58 

85.4 

5264.81 

3.22 

-3.13 

-32.01 

-4.48 

-4.48 

49.2 

5284.11 

2.88 

-3.11 

-32.04 

-4.57 

-4.57 

60.3 

5325.55 

3.21 

-3.16 

-32.03 

-4.55 

-4.55 

45.2 

5414.07 

3.21 

-3.58 

-32.02 

-4.52 

-4.52 

29.8 

5425.26 

3.20 

-3.22 

-32.04 

-4.56 

-4.56 

43.4 

5991.38 

3.15 

-3.55 

-32.05 

-4.58 

-4.58 

32.1 

6239.95 

3.89 

-3.46 

-32.00 

-4.55 

-4.55 

13.2 

6247.56 

3.89 

-2.30 

-32.00 

-4.58 

-4.58 

55.0 

6369.46 

2.89 

-4.11 

-32.06 

-4.59 

-4.59 

20.7 

6432.68 

2.89 

-3.57 

-32.07 

-4.56 

-4.56 

43.2 

6456.38 

3.90 

-2.05 

-32.00 

-4.59 

-4.59 

66.9 


Note. — q 

'’-values are fromIBlackwell et al.l 11979lll982allbll: 

O’Brian et al. 

Il991fl. andfBard et al.| Il99ll) for Fe I and from 

Melendez Barbuyl 120091) for Fe II. 
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TABLE 5 

Final stellar parameters. 


HD, 

BD 

Ted 

K 

log g 

[Fe/H] 

km 

mass 

Mq 

Fe I - Fe II 
LTE 

Fe I - Fe II 
NLTE 

Ypel 

Ypell 

19373 

6045 ± 80 

4.24 ±0.05 

Benchmark stars 

0.10 ±0.05 1.2 1.11 

0.06 ±0.09 

0.06 ±0.09 

26 

15 

22484 

6000 ±100 

4.07 ±0.05 

0.01 ±0.04 

1.1 

1.18 

-0.02 ±0.07 

0.00 ±0.07 

27 

16 

22879 

5800 i 90 

4.29 ±0.07 

-0.84 ±0.07 

1.0 

0.75 

-0.06 ±0.08 

-0.03 ±0.08 

23 

14 

30562 

5900 ± 85 

4.08 ±0.05 

0.17 ±0.08 

1.3 

1.12 

0.02 ±0.09 

0.02 ±0.09 

26 

16 

34411 

5850 ±100 

4.23 ±0.05 

0.01 ±0.03 

1.2 

1.10 

-0.01 ±0.05 

-0.02 ±0.05 

26 

15 

49933 

6600 ± 80 

4.15 ±0.05 

-0.47 ±0.07 

1.7 

1.11 

-0.05 ±0.08 

-0.04 ±0.08 

26 

15 

59374 

5850 ± 55 

4.38 ±0.09 

-0.88 ±0.05 

1.2 

0.75 

-0.06 ±0.09 

-0.02 ±0.09 

24 

15 

59984 

5930 ±100 

4.02 ±0.06 

-0.69 ±0.07 

1.4 

0.89 

-0.06 ±0.10 

-0.06 ±0.10 

23 

18 

64090 

5400 ± 70 

4.70 ±0.08 

-1.73 ±0.07 

0.7 

0.62 

-0.02 ±0.09 

-0.02 ±0.09 

22 

11 

69897 

6240 ± 70 

4.24 ±0.05 

-0.25 ±0.04 

1.4 

1.15 

-0.06 ±0.07 

-0.04 ±0.07 

28 

15 

84937 

6350 ± 85 

4.09 ±0.08 

-2.12 ±0.07 

1.7 

0.76 

-0.06 ±0.11 

0.00 ±0.12 

12 

7 

94028 

5970 ±130 

4.33 ±0.08 

-1.47 ±0.04 

1.3 

0.70 

-0.06 ±0.07 

-0.04 ±0.07 

20 

15 

102870 

6170 ± 80 

4.14 ±0.04 

0.11 ±0.06 

1.5 

1.35 

-0.03 ±0.07 

-0.01 ±0.07 

23 

14 

103095 

5130 ± 65 

4.66 ±0.08 

-1.26 ±0.08 

0.9 

0.60 

0.01 ±0.11 

0.01 ±0.11 

22 

8 

105755 

5800 zb 55 

4.05 ±0.09 

-0.73 ±0.05 

1.2 

0.85 

-0.01 ±0.06 

0.00 ±0.06 

23 

15 

114710 

6090 ± 80 

4.47 ±0.05 

0.06 ±0.06 

1.1 

1.17 

0.03 ±0.07 

0.04 ±0.06 

26 

15 

134169 

5890 ± 80 

4.02 ±0.07 

-0.78 ±0.07 

1.2 

0.87 

0.01 ±0.09 

0.06 ±0.09 

26 

18 

140283 

5780 ± 55 

3.70 ±0.07 

-2.46 ±0.07 

1.6 

0.80 

-0.06 ±0.09 

-0.02 ±0.09 

20 

14 

142091 

4810 ± 65 

3.12 ±0.06 

-0.07 ±0.10 

1.2 

1.16 

0.05 ±0.13 

0.05 ±0.12 

20 

13 

-1-66° 0268 

5300 ± 80 

4.72 ±0.11 

-2.06 ±0.15 

0.6 

0.60 

-0.03 ±0.16 

-0.03 ±0.16 

21 

9 

24289 

5980 

Stars with spectroscopic parameters 

3.71 -1.94 ±0.17 1.1 0.83 -0.10 ±0.19 

-0.03 ±0.19 

16 

10 

30743 

6450 

4.20 

-0.44 ±0.07 

1.8 

1.08 

-0.03 ±0.09 

0.01 ±0.09 

26 

18 

43318 

6250 

3.92 

-0.19 ±0.08 

1.7 

1.26 

-0.05 ±0.13 

-0.01 ±0.09 

28 

15 

45067 

5960 

3.94 

-0.16 ±0.06 

1.5 

1.18 

-0.03 ±0.07 

-0.01 ±0.07 

24 

16 

45205 

5790 

4.08 

-0.87 ±0.03 

1.1 

0.82 

-0.08 ±0.07 

-0.05±0.07 

23 

16 

52711 

5900 

4.33 

-0.21 ±0.05 

1.2 

0.97 

0.03 ±0.07 

0.05 ±0.07 

26 

16 

58855 

6410 

4.32 

-0.29 ±0.05 

1.6 

1.15 

-0.01 ±0.08 

0.01 ±0.07 

26 

15 

62301 

5840 

4.09 

-0.70 ±0.04 

1.3 

0.85 

-0.07 ±0.08 

-0.04 ±0.08 

29 

16 

74000 

6225 

4.13 

-1.97 ±0.07 

1.3 

0.76 

-0.08 ±0.10 

-0.02 ±0.10 

15 

7 

76932 

5870 

4.10 

-0.98 ±0.05 

1.3 


-0.01 ±0.07 

0.03 ±0.08 

27 

16 

82943 

5970 

4.37 

0.19 ±0.04 

1.2 

1.17 

0.01 ±0.06 

0.01 ±0.06 

25 

13 

89744 

6280 

3.97 

0.13 ±0.03 

1.7 

1.50 

0.01 ±0.05 

0.02 ±0.06 

26 

13 

90839 

6195 

4.38 

-0.18 ±0.05 

1.4 

1.10 

0.04 ±0.07 

0.06 ±0.07 

28 

17 

92855 

6020 

4.36 

-0.12 ±0.03 

1.3 

1.08 

0.00 ±0.06 

0.02 ±0.05 

24 

11 

99984 

6190 

3.72 

-0.38 ±0.04 

1.8 

1.33 

-0.01 ±0.07 

0.03 ±0.07 

24 

15 

100563 

6460 

4.32 

0.06 ±0.08 

1.6 

1.30 

0.00 ±0.10 

0.02 ±0.10 

23 

10 

106516 

6300 

4.44 

-0.73 ±0.06 

1.5 

0.95 

-0.03 ±0.08 

0.00 ±0.08 

22 

15 

108177 

6100 

4.22 

-1.67±0.05 

1.1 

0.72 

-0.08 ±0.13 

-0.06 ±0.14 

16 

5 

110897 

5920 

4.41 

-0.57 ±0.04 

1.2 

0.85 

0.03 ±0.05 

0.05 ±0.05 

29 

18 

115617 

5490 

4.40 

-0.10 ±0.05 

1.1 

0.93 

-0.04 ±0.08 

-0.05 ±0.08 

26 

13 

134088 

5730 

4.46 

-0.80 ±0.05 

1.1 

0.75 

-0.02 ±0.07 

0.00 ±0.07 

25 

15 

138776 

5650 

4.30 

0.24 ±0.05 

1.3 

0.93 

0.02 ±0.10 

0.00 ±0.10 

21 

11 

142373 

5830 

3.96 

-0.54 ±0.05 

1.4 

0.95 

-0.02 ±0.07 

-0.02 ±0.07 

28 

17 

-4° 3208 

6390 

4.08 

-2.20 ±0.09 

1.4 

0.77 

-0.08 ±0.11 

0.00 ±0.11 

16 

11 

-13° 3442 

6400 

3.95 

-2.62 ±0.09 

1.4 

0.79 

-0.14 ±0.11 

0.00 ±0.11 

8 

5 

-1-7° 4841 

6130 

4.15 

-1.46 ±0.05 

1.3 

0.79 

0.00 ±0.09 

0.02 ±0.08 

21 

17 

-1-9° 0352 

6150 

4.25 

-2.09 ±0.04 

1.3 

0.70 

-0.01 ±0.07 

0.03 ±0.07 

15 

6 

-1-24° 1676 

6210 

3.90 

-2.44 ±0.09 

1.5 

0.78 

-0.06 ±0.12 

0.04 ±0.12 

12 

5 

-1-29° 2091 

5860 

4.67 

-1.91 ±0.08 

0.8 

0.70 

-0.05 ±0.10 

-0.05 ±0.10 

18 

10 

-1-37° 1458 

5500 

3.70 

-1.95 ±0.09 

1.0 


-0.07 ±0.11 

-0.05 ±0.11 

21 

14 

G090-003 

6010 

3.90 

-2.04 ±0.06 

1.3 

0.78 

-0.06 ±0.10 

-0.02 ±0.10 

17 
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